Collapse models with non-white noises II: particle-density coupled noises 
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^ . We continue the analysis of models of spontaneous wave function collapse with stochastic 

^ ' dynamics driven by non-white Gaussian noise. We specialize to a model in which a classical 

^vg I "noise" field, with specified autocorrelator, is coupled to a local nonrelativistic particle den- 

,_^ , sity. We derive general results in this model for the rates of density matrix diagonalization 

O,. and of state vector reduction, and show that (in the absence of decoherence) both processes 

■y ' are governed by essentially the same rate parameters. As an alternative route to our re- 

Cd ' duction results, we also derive the Fokker-Planck equations that correspond to the initial 

^^, stochastic Schrodinger equation. For specific models of the noise autocorrelator, including 
ones motivated by the structure of thermal Green's functions, we discuss the qualitative and 

^ ' qantitative dependence on model parameters, with particular emphasis on possible cosmo- 

._!. ' logical sources of the noise field. 
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In an earlier paper [l[, hereafter referred to as (I), we presented a detailed analysis of stochastic 



H ' models for state vector collapse driven by Gaussian non-white noise. In particular, we showed 

that a perturbation expansion in the noise strength parameter yj^ permits the explicit calculation 
of consequences of the model, in parallel with standard results obtained by the Ito calculus in 
the white noise case. In (I) the noise couplings were introduced in generic form, subject to the 
assumption that the noise correlator has a positive definite structure in the large time limit. As 
we shall see, this positivity assumption is overly restrictive, and does not apply to the physically 
interesting case of thermal noise, where the spatial Fourier transform of the noise correlator is 
oscillatory in time. Our aim in this paper is to specialize the discussion of (I) to the physically 
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interesting case of a particle density-coupled classical noise field, and then within the context of 
this model, to give a generalized analysis of density matrix diagonalization, state vector reduction, 
and constraints on model parameters. We then turn to the question of whether the noise field 
postulated in stochastic reduction models can be realized as a cosmological field. (For a section- 
by-section brief summary of the contents of this paper, the reader should turn to the Summary 
and Discussion given in Sec. VII.) 

Our starting point in (I) was a diffusion process for the wave function in Hilbert space having 
the form (with h = 1, with the constant complex coupling factor ^ introduced in (I) set equal to 
1, and with the state vector denoted here by {tp)), 

N 

-iH + ^^A,u;i(t) + O 



dm)) 



^^ .„ , , , ^-_.,^, , ^ m)) . (1) 

Here H is the standard quantum Hamiltonian of the system, Ai are commuting self-adjoint oper- 
ators, 7 is a positive coupling constant, and O is a linear operator yet to be defined. The noises 
Wi{t) are real Gaussian random processes, whose mean and correlation functions are, respectively 

EK(t)]=0, E[w^ih)wj{t2)]=Dijih,t2) . (2) 

We will now specialize the discussion to the case in which the index i is the spatial coordinate 
X, and the operator Ai is a particle density M{x), which, for a many body system composed of 
distinguishable particles with couplings rrii and coordinate operators iji, is given by 

M{x)=Y,m6Hx-qi) . (3) 

i 

(We have chosen a notation appropriate to the case in which the density M is a mass density, 
but ([3]) also describes other forms of coupling to particle densities, such as to the baryon number, 
lepton number, or isospin densities, with m,j the appropriate coupling constants.) An important 
property of the density operator of ([3]) is that when integrated over space it reduces to a c-number 
that commutes with all operators, 

/ (fxM{x) = / (fx'^mi6^{x- qi) = ^m-j. (4) 

i i 

Hence a noise coupling to the density operator can be permitted to have a nonzero expectation, 
since this will only contribute a constant term to the effective Hamiltonian on the right of ([1]). So 
we will assume that, corresponding to (l3|), the noises Wi{t) of (I) form a classical noise field, which 
we shall denote by (j){x,t), with mean and autocorrelation 

E[(t){x,t)] = cf>o, E[{cP{x,ti)-cf>o)){cf>iy,t2)-(t>o)] = D{x-y,ti-t2) . (5) 



Here, in assuming a constant expectation (f)Q and in writing the arguments of D, we have built 
in an assumption of space and time translation invariance; we shall also assume spatial inversion 
invariance, so that D{x,t) = D{—x,t). Thus, with this specialization of the noise structure of (I), 
the diffusion process in Hilbert space of ([1]) becomes 

dm)) 



dt 



-iH + ^ I d^xM{x)<i){x,t) + O 



Hit)) . (6) 



In most of what follows, we will neglect the Hamiltonian term in ([6]), focusing on effects that arise 
from the action of the stochastic term. 

Because it uses real-valued noise, ^ does not preserve the norm of the wave function, and 
this is where the operator O enters. In (I), through detailed calculations that we shall not repeat, 
we show that O is fixed by the requirements of (i) state vector normalization, and (ii) a linear 
evolution equation for the density matrix 

p{t) = Emt))m)\] , (7) 

which guarantees that superluminal signaling cannot occur. Relation ([7]) guarantees also the posi- 
tivity of p{t) throughout time. Determining the structure of O leads to three equations from (I), 
which are exact to order 7, and which when specialized to the case of a density coupled noise, form 
the starting point for our analysis here. 

The first of the needed equations describes the density matrix time evolution, as given by (53) 
of (I), 

H r^ r r 



^-p{t) = -i[H, p{t)]+j J ds J d-'x J d-'y[M{x) p{t) M{y, s - t) + M{y, s - t) p{t) M{x) 
-M{x)M{y, s - t)p{t) - p{t)M{y, s - t)M{x)]D{x - y,t - s) , 

(8) 

with [see (47) of (I)] 

M{y, s-t) = e*^(^-*)M(y )e-^^(^"*) . (9) 

When H = this simplifies to read [see (19) of (I)] 

j^pit) =7 / d^x f d^y[M{x) p{t) M{y) + M{y) p{t) M{x) 



-M{x)M{y)p{t) - p{t)M{y)M{x)]F{x - y, t) 



(10) 



where we have defined 

F{x-y,t) = / dsD{x-y,t- s). (11) 

Jo 

To state the second equation, let us define the expectation {0)t = {ij{t)\0\ip{t)) for any operator 

O. Then when H = the time evolution of the stochastic expectation of the variance V^(t) = 

{A'^)t — {A)f of any operator A that commutes with the mass density for all x, given by (23) and 

(24) of (I), becomes 

j^nVAit)] = - 87 I d'x I d'yE[{{Mix)-{M{x))t)A)t{{M{y)-{M{y))t)A)t]Fix-y,t). (12) 

The final equation that we need describes the time evolution of the state vector |'0(O)) as specified 
in (40), (51), and (52) of (I), which combined become 

dm)) 



dt 
with the self-adjoint operator B given by 



iH + ^ j d\[M{x)-{M{x))Mx,t) + j{B-{B)t) 



mt)) , (13) 



B = -2j d^xj d^yF{x - y, t) [M(f ) - (M(x))i] [M(y) - {Miy))t] . (14) 

The alternative form of this equation given in (35) and (37) of (I) differs only by a change of 
measure for the noise, and makes identical physical predictions. 

As is easily checked, an important consequence of the fact that the spatial integral of M(x) is 
a c- number (c.f. (jj])) is that the noise field expectation (J)q makes no contribution to the order ^/7 
term in (fT3|) . and that a space-independent constant in F{x — y,t) makes no contribution to P^ . 
p^ . and (dH). That is, we can replace F{x — y,t) by the subtracted function 

F{x-y,t)-m , (15) 

for an arbitrary function £,{t), with no effect on the equations; only the nonzero spatial Fourier 
components of F(x — y,t) are significant for our analysis. In particular, for S,{t) = F{0,t), this 
invariance implies that we are free to replace F{x — y,t) by the subtracted function F{x — y,t) — 
F{0,t), which has a spatial Fourier transform with improved convergence at small wave numbers. 
There has recently been a spirited debate [4] over whether stochastic reduction models can 
be made relativistically invariant. We remark in this context that the noise coupling of ([1]) can 
be obtained in a number of ways as the non-relativistic limit of relativistically invariant, anti- 
self-adjoint coupling actions involving scalar, vector, or tensor fields. (An anti-self-adjoint action 
is required to give a real noise term in the Schrodinger equation; we will not attempt here a 



fundamental justification of this phenomenologically- motivated choice of Hermiticity structure.) 
When the noise coupling is introduced as the nonrelativistic limit of a relativistic action, relativistic 
invariance of the stochastic reduction model is broken not by the noise coupling, but by the assumed 
autocorrelator of the noise field (j){x,t). For example, if the noise field has a cosmological origin, 
its autocorrelator might be expected to refer preferentially to either the Lorentz frame in which 
the cosmological background radiation is isotropic, or to the galactic rest frame. A topic for future 
work will be to investigate whether an effective anti-self-adjoint coupling action can arise naturally 
in a non-equilibrium cosmology, or requires an explicitly non- unitary pre-quantum dynamics. 

II. DENSITY MATRIX DIAGONALIZATION 

We begin our analysis by considering the consequences of ()10p for coordinate off-diagonal ma- 
trix elements of the density matrix, when the Hamiltonian evolution is neglected. Taking the 
matrix element of ([TO]) between states ({r|}| and |{r|}), we get a differential equation for the time 
dependence of the matrix element of p, which can be immediately integrated to give 

{{r}}\p{tWh) = e-^'^'HiriMpiOMrn) , (16) 

with the integrated rate T{t) given by 

r(t) = j (fix (fiy / dsF{x-y,s)[mi{x) - m2{x)][mi{y) - m2{y)] , (17) 
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where mi^2 are the eigenvalues of the operator M(x) when acting on the respective states \{ff, ' }), 





iTii{x) = y^mi6'^{x — fl) , m2{x) = y^mi6'^{x — f^) . (18) 



Substituting (J18l) into (J17p and carrying out the x and y integrals using the delta functions, we 
obtain 

•^0 i j 

We now review a number of useful features of this formula (many of which, in a slightly different 
notation, are familiar from the stochastic reduction literature). First of all, as already pointed out 
in Sec. 1, (|19|) is unchanged when we replace F{f, s) by the subtracted function F{f, s) — F{0, s). 
Secondly, suppose that rj=rf = f/ for some particle with index /. Then the contribution of this 
particle to the double sum in p^ is 

2m?[F(0, s) -F(0, s)] +2m/ ^ mj[F{ri-rj, s) +F{fi-ff, s)-F{fi-ff, s) -F(rf -r/, s)] = . 

(20) 



So only particles that have different coordinates in the groups 1 and 2 contribute to the sum. 

Third, suppose that for large separations r, relative to some correlation scale re, the function 
F(r',s) asymptotically approaches a constant (which can be zero or nonzero). Then if there are 
two particles /, J such that fj' — fy are all large enough relative to re to be in the asymptotic 
regime for F, the cross terms in the double sum linking these two particles do not contribute. This 
means that if the particles form a set of K widely spaced bunches on the scale of re, with the 
particles of group 2 displaced with respect to those of group 1 by distances of order re, the formula 
for r(i) splits into a sum 

T{t) = J2^^{t) , (21) 

fc=i 

with T^{t) computed entirely within the A;th bunch. 

Fourth, let us take group 1 to be a collection of particles that are very closely spaced on the 
scale of re, and suppose that the particles of group 2 are all displaced by a common vector R with 
respect to those of group 1. In this case, T{t) is approximated by the formula 

T{t) =. 27 / ds\ymi\ [F{6, s) - F{R, s)] , (22) 



' ds I y^ rui 

•^0 V ^ / 



which is the formula that would be obtained if there were only one particle of mass ^^ rrii at 
the center of mass of the group. The above formulas display the amplification mechanism typical 
of collapse models: when particles interact to form a macro-object, the collapses on the single 
particles add up in such a way that the center of mass of the object collapses each time a single 
particle does. This is the reason why these models can account both for the quantum properties 
of microscopic systems and for the classical properties of macroscopic objects. 

Fifth, let us again take group 1 to be a collection of particles that are very closely spaced on 
the scale of re, but now suppose that the particles of group 2 are displaced by random amounts, 
with an average magnitude of displacement R with respect to those of group 1. When the function 
F{f,s) only depends on the magnitude \r\ of the displacement vector, so that F{f,s) = F[\r\,s] 
then r(t) is approximated by the formula 



rit)^jJ^dslS2mA [FiO,s) + {{Fi\f^-f^\,s)))N-2{{Fi\fl-r^ls)))N] , (23) 

where {{■■■))n denotes the average over the ensemble of particles; when R > re (|23p is further 
approximated by 

r(t):^7^*dJj]m,j [F{0,s)-F{R,s)] , (24) 



which is one half of the T(t) given by the center of mass formula (j22p for the corresponding 
magnitude of R. 

Finally, in many cases of interest F{x — y, s) can be written as a sum or integral over factors 
referring to x and y separately, 



F{x-y,s) = ^J^{a,s)f{a,x)f{a,y) , 



(25) 



with J^ an appropriate weighting function, and a a shorthand for any combination of discrete and 
continuous variables. Substitution of (f25]l into (flTl) gives 

rt 
/O 
with 



r(t) = 7 / dsY,Hc^,s)g{a,{fl},{ri}f 
Jo 



(26) 



g{a,{rl},{f^}) = I d^xf{a,x)[mi{x) - m2{x)] =Y,^^[f{c^,rl) - f{a,r^)] . (27) 



III. STATE VECTOR REDUCTION 



We proceed next to apply (|T2|) to the problem of reduction of a state vector constructed as a 
superposition of position eigenstates |{r^}) , j = 1,...,N. Our first step is to rewrite ([12]) in a 
more useful form by setting A = B + C, with B and C operators that commute with each other 
and with the mass density, and subtracting off (J12p as written for B and C alone, which gives 

d_ 
di' 



E[{BC)t-{B)t{C)t] = - 87 I d'x I d'yE[{{M{x)-{M{x))t)B)t{iMiy)-{M{y))t)C)t]F{x-y,t) 

(28) 



Using the fact that B and C can be arbitrary operator functions of the particle coordinate operators, 
by making the specific choices -B = ]^^ 6^{ui — qi) and C = J^^ Si^{wi — iji), we get 



dt 



n sHm - n,) m{w,})\' - m{uimH{^i})\' 



87 / d-'xd-'yE 



\i,i{u,})\'m{we})\' 



Y,mj[6^{x-Uj) - \i;j{x)\^]J2mk[6^{y-Wk) - \My)(' 



F{x-y,t) , 



(29) 



where we have introduced the definition 



(30) 

Let us now speciahze (|29]) to the case of a wave function which is the superposition of N distinct 
localized groups of particles, by writing 

N 

i^m]) = {{m^{t))=Y,c^jj{5\z,-rif/' , 

j=i I 

N 

j=i e. 

(31) 

with pj = |ajp and with normalization of the wave function implying that "^jPj = 1- (By 
the square root of a delta function, we mean a Gaussian wave packet which is sharply localized, 
with a modulus squared that integrates to unity.) Substituting ([3T]) into ([29]) . and integrating in 
d{w£} around {r^ } and in d{'Ui} around {r* }, we get an equation for the time evolution of the 
occupation probabilities pj of the corresponding states Yie ^^{^e ~ rjY'"^ with label J that appear 
in the superposition, 
d 



-E[6mlPai -PmPl] = - 87E 



j k [ R,S 



PMPL }^ y^mjmk{ F(f/^ - r;:\t) + }^PRPsF{rf - r^ ,t) 



R S 

(32) 
Specializing this further to the two group case with N = 2, taking M = L = 1 and doing some 
algebraic rearrangement using the fact that the sum of the probabilities is pi +p2 = Ij we get 

^EbiP2] = - 8jE[pIpI] Y E ^j^k{F{fl - rlt) + F{rj - ri,t) 

j k 

-F{rj-rlt)-F{rj-rlt)} . 

(33) 

We can now use ([33|) to derive upper and lower bounds to the reduction rate, as follows. To 
obtain an upper bound, we use the inequality 

¥.\plp\]>W.\p^P2f , (34) 



and the assumption that the integrand of T{t) in p9]) is positive for all s, to rewrite ([33|) as 
J^nPiP2] < - 8jE[p,p2? Yl E ^J^k{F{rl - rlt) + F(f/ - rlt) 

j k 



-F{rl-flt)-F{r^-rlt)] , 



(35) 



giving a differential inequality that can be integrated to give an upper bound on the reduction rate 

Eb.Wp.Wl£»2M , (36, 

To get a lower bound, we note that since the probabilities pi and p2 obey pi + P2 = 1, we have 
PiP2 = Pi(l — Pi) < 1/4, and so 

E[py2] < E[piP2]/4 . (37) 



Again assuming that the integrand of (J19p is positive for all s, this gives a differential inequality 
that can be integrated to give the lower bound 

E[Pi(i)P2(i)]>lE[m(0)p2(0)]exp[-2r(t)] . (38) 

Thus we see that in our model of a Schrodinger equation modified solely by a real noise process, the 
upper and lower bounds on the reduction rate involve (under the uniform positivity assumption) the 
same integrated rate function T{t) as appears in the decay of the off-diagonal density matrix element 
({^/}|p(^)|{^/})- Of course, in realistic applications, the rate for density matrix diagonalization is 
expected to receive much larger contributions from decoherence processes, which can be modeled by 
imaginary noise terms in the Schrodinger equation that do not contribute to state vector reduction. 
Although the upper and lower bounds are governed by the same integrated rate function, they have 
very different functional dependencies: the upper bound depends on the inverse of Tit), whereas 
the lower bound is a negative exponential in T{t). Solvable models [3], [J|, and Appendix D, show 
that in fact the actual decay of the variance is exponential, rather than power law, indicating that 
the lower bound of ([38]) gives the better estimate.^ 



^ In the example calculated in Appendix D, the actual variance decay is --^ e ^^*' . A simple example shows how 
an exponential decay of the variance can agree with the inequalities used to get the upper and lower bounds. 
If E[pi(i)p2(t)] ^ E[pi(0)p2(0)]exp(-r(t)), then (d/dt)E[pi(i)p2(t)] = -r'(i)E[pi(i)p2(i)], whereas ^ implies 
that (d/rff)E[pi(i)p2(t)] = -9>T' {t)'E[pl{t)pl{t)], and so we must have E[pi(i)p(t)] = SEfpfpj]- Suppose now that 
pi{t)p2{t) = with probability 1 — exp(— r(i)), and pi{t)p2{t) = 1/8 with probability exp(— r(i)). We then have 
E[pi{t)p2{t)] = 8E[pi{t)pl{t)] = (l/8)exp(-r(t)). However, E[pi(t)p2{t)f = (1/64) exp(-2r(t)), which for large 
r(f) is much smaller than E[pf{t)p2{t)] = (1/64) exp(— r(f)), and so the inequality of (|34p . which was used to get 
the upper bound, is far from being saturated, while by construction, the inequality of (|37|l . which was used to get 
the lower bound, is saturated. 



10 



In the general N group case, although we have not derived rigorous bounds, we can get estimates 
similar to the two-group case by setting Af = L in (I32p . giving 



-E[pl{1 - Pl)] = - 8jE 



j k [ R,S 



ri,t) 



R S 



(39) 



Suppose now that the stochastic process brings the probabilities close to a corner of their domain, 
where for some M ^ L the probability pM is close to unity, and thus all the other probabilities are 
small. The right-hand side of (j39j) then contains terms of second degree in small quantities, given 
by selecting the terms with R = M and S = M in the sums, plus remaining terms that are third 
degree in small quantities. The second degree terms contribute 

^E\pl] c^ - 8jE[pI] Y^ Y1 ^o^k\F{rj^ - rt . t) + F{f^^^ - ff , t) 
j k I 

-F{r^-rt.t)-F{r^-rt\t)\ , (40) 

which has a structure similar to ([33ll for the two-group case. Using the inequality 

E\pl]>E\pLf , (41) 

defining V^^{t) by 

r^"'W =7 r^s j;^m,mfc[F(rf -rt,.) + F(r-;^-ff ,.)-F(rf -ff ,.)-F(r;*^-rt,5^ , 
•^° j k 

(42) 



and assuming positivity of the integrand of (j42p . we get a differential inequality that can be 
integrated to give an upper bound on the decay rate, 

.,.„„_||^ . ,3, 

Similarly, from (j40p we can also get a lower bound on the decay rate, 

E[pi(t)]>E[pL(0)]exp[-2r^*^(t)] . (44) 

Thus, near the corner where pM — li the other pi decay to zero, with the slowest rate of decrease 
corresponding to the smallest value of r for L ^ M. 
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To conclude this section, we note that when F{x — y,t) has the factorized form given in ()25p . 
then (I32p takes the form 



-t-K[5mlPm - PmPl] = - 87E 



PMPL 



^J^{a,t)^m,[fia,fl^)-Y,PRfi(^,rP)] 



R 



Y^ mk[f{a, rt^) - ^psfia, ff )] 



while (l42l) becomes 



■^0 a i 

=^ f dsY,H»,t)9ia,{r[},{rnf , 
-'0 



and T(t) = T (t) is the specialization of this formula to L = 1, M = 2. 



(45) 



(46) 



IV. FOKKER-PLANCK EQUATION 

As a complement to the methods used in the preceding sections, we derive the Fokker-Planck 
equation for the non-white noise model, and use it to rederive ()32p . We again restrict ourselves 
to the case when the Hamiltonian H is zero, which allows all equations to be diagonalized in 
coordinate representation. Starting from ([T3D and substituting 

N 
L=l 

with |{r/'}) sharply localized wave-packet states [c.f. ([3T]) ]. we find that the coefficients a^ obey 
the equation of motion 

d 



dt 



ai = chlXl , 



(48) 



12 



with Xl given by 



Xl=ViY. 



mi[(t){rt^ 



t)-Y,PR^ir,'',t)] 



R 



-7 



Y^Y.'^m. 



« J 



Firt-rj- 



-2Y,PrF{ 



J 



t) + 2Y,PRPsF{r^-r^,t) 
R,S 

,t)-Y,PRHr^''-fj',t) 



R 



R 



Since Xj^ is real, and pL = a*^aL, we correspondingly have 

d 

—PL = 2plXl . 
at 



(49) 



(50) 



In order to derive the Fokker-Planck equation, we have to evaluate E[(d/(it)/({p/j})] for an 
arbitrary function / of the set of probabilities {p_r}, keeping terms through order 7. On using the 
chain rule we have 



E 



f{{PR}) 



E 



E 



a/({TO}) d 

dps dt 



PS 



(51) 



On substituting (j50p and (j49p for (d/dt)ps, we encounter two types of terms. Terms of the form 

'dfiiPR}) 



^Y^E 



dps 



-A, 



(52) 



can be read off directly from the term proportional to —7 in Xs, while terms of type 

'dfiiPR}) 



VtE^ 



i,S,L 



dps 



-Bs^L(t>{r|',t 



(53) 



are evaluated using the Furutsu-Novikov formula, which approximated by using pr{s) = PR{t) + 
0(^7), takes the form 



E 



dfiiPR}) 
dps 



Bs^LHfl',t) 



j,S,T 



d 



dpT y dps 



BsiL 



dpi 



d4>{ff , t) 



(54) 



The needed derivative of pT can be read off directly from the ^77 term in Xt', substituting this, 
and doing much algebra, one finds that all first derivatives of / with respect to the probabilities 
cancel exactly, leaving finally the compact expression 



d 
di 



nf{{pR})] = 47 Y. ^^""^ - ^/' *)"^*"^jE 

M,T,i,i,Q,S 



d^f{{PR}) 
dpNdpM 



PTPMihiQ - Pq){StS - Ps) 



(55) 



d_ 

IE 
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Introducing the probability density P({p/j},t), which includes as a factor the constraint 
5{Y1,lPl — 1) requiring that the probabilities sum to unity, we can also write the expectation 
of {d/dt)f{{pR}) as 

nf{{VB})] = \{jdpL^-^^^^f{{VR}) . (56) 

Li 

Comparison of this expression with (j55p . as rearranged by two integrations by parts (the surface 
terms when any probability is or 1 do not to contribute; see below), one gets the Fokker-Planck 
equation, 

^^=Ea^[-^«T(te))^(te).')l . (57) 

with 

Amt =4'ypMPT ^ F{fl^ - rf , t)mimj {5mq - Pq){Sts - Ps) 

i,j,Q,S 

=^1PMPT Y. rn,mjPQPs[F{rf - rj, t) + F{r^ - r^ , t) - F{r^ - rj, t) - F(r;^ - ff , t) 
i,j,Q,S 

(58) 



This equation is a specific case of a general Fokker-Planck equation written down by Pearle [5[ as 
the basis for a general class of objective reduction models. When F{x — y, t) takes the factorized 
form of (|25p . Amt can be rewritten as 



Amt = '^IPMPT y^ -^(a, t)(j)M4'T , (59) 

a 

with 

cl>M = Y.miPQ[f{a,rt')-fic^,r^)] ■ (60) 

i,Q 

We see that in addition to vanishing when either pM = or pT = 0, Amt vanishes when either 
PM = 1 or PT = 1, because (pM vanishes when pM = 1- This is why the integrations by parts leading 
to the Fokker-Planck equation produce no surface terms, and also why the Fokker-Planck equation 
of (|57|) satisfies the criteria that Pearle [5| formulated for getting a Fokker-Planck equation that 
leads to state vector reduction with Born rule probabilities. 

As an application of ([55]) . if we substitute /{{pr}) = Pl we find, since d^pL/{dpNdpM) = 0, 
that 

Jt^iPL] = . (61) 



14 



Similarly, if we substitute /{{pr}) = SklPl — PkPl-, we find, using 

52 



-[^KLPL -PKPh] = -[SmK^NL + SmL^Nk] , (62) 



dpAidpN 

that (I55p yields (I32p . More generally, (I55p and the corresponding Fokker-Planck equation of ()57p 
allow one to calculate the time evolution of a general function /({p_r}) of the probabilities. 

V. NOISE EFFECTS: ENERGY PRODUCTION AND RADIATION BY ATOMS 

The noise coupling postulated in Sec. 1 as the origin of state vector reduction has other physical 
effects, that serve to place upper bounds on the noise coupling strength 7. We focus in this section 
in particular on energy production, and gamma radiation from atoms, which place particularly 
stringent bounds on the model parameters. 

A. Energy production 

To calculate the mean rate of energy production, we have to evaluate {d/dt)TTHp(t) = 
TrH{d/dt)p{t). From ^ we find, by repeated cyclic permutation under the trace, that 

jTTHp{t) = -^Jd^xld^yl dsD{x-y,t-s)TT{[[H,M{x)],M{y,s-tMt)) . (63) 

This equation is exact through order 7. We now make the Markovian approximation, of assuming 
that we can ignore the "memory effect" associated with the characteristic decay time of the noise 
correlator D{x — y,t — s), by replacing M(y,s — t) by M{y,0) = M{y). For white noise, where 
D{x — y,t — s) = G{x — y)5{t — s), the Markovian approximation is exact; for non- white thermal 
noises, it should be a good approximation when the energy at the peak of the noise spectrum 
is much higher than the typical kinetic energies of the particles to which the noise couples (see 
Appendix A). With this approximation, ()63p simplifies to 

j^l)cHp{t) = --fJd^xld^yF{x-y,t)Tt{[[H,M{x)],M{yMt)) , (64) 



where we have made use of the definition (jlip . 

Let us now assume that H is the nonrelativistic Hamiltonian for a collection of particles inter- 
acting through a general velocity-independent potential, 



^ = E^ + ^(i^a) , (65) 
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while M(x) has the form of ^ (in the mass-coupled noise case, rrii = Mi) and F{x — y,t) has the 
factor decomposition of ()25p . Then carrying out the x and y integrals, (I64p becomes 



d 



-TtHp{t) = -^^J^{a,t)^mim,Tt{[[H, f{a,(u)], f{a,qjMt)) 



dt 



(66) 



«.j 



The commutators appearing in (j66p are easily evaluated, 

r2 



[HJia,^ 



P 



E^ + ^(i^^»'/(«'9^) 



Pi 

2M, 
—i 



■,f{a,qi) 



2Mi 



\Pi ■ ^qrfioi: Qi) + Vg;/(a, gj) • Pi] 



giving 



[[H,f{a,ql)],f{a,q,)] = -6,, — [V,^J{a,q^f 



Substituting this into (j66p . we obtain finally 
d 



dt 



Tri/p(t) = 7 5] ^(a, t) 5^ TFTr ([V,-;J(a, gl)] Vt)) 



(67) 



(68) 



(69) 



A further simplification of this result can be achieved by using the Fourier transform represen- 
tation of F(x — y,t), which (recalling the assumed spatial inversion invariance) takes the form 

d^k 



F{x-y,t) 



(2^ 
_fk_ 

(2^ 



cos fc- {x-y))F{k,t) 



cos(/c • x) cos{k ■ y) + sin(A: • x) sin(A; • y)\F{k, t) 



(70) 



This has the general structure of ([25|) . with ^^ corresponding to j d^k{2TT)^'^Y^^^]^, with n a 
discrete index distinguishing between the sine and cosine modes, that is, J-{a, t) = F{k, t) for both 
n = 0,1, and f{k, n = 0,x) = cos{k • x) and f{k, n = l,x) = sin{k • x). Substituting ([70]) into ([69|) 
then gives 



'^TrHp{t)=j[-^Fik,t)Y,^Tr 



dt 



(27r)3 

d^k r 

(2^ 



fc^f cos^(A; • (ji) -|-sin^(A; • qi))p{t) 



pFik,t)Y, 



M, 



(71) 
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where in the final step we have used Trp{t) = 1. Thus, in the Markovian approximation, we get a 
simple formula for the energy production rate, expressed entirely in terms of the Fourier transform 
of F(x — y,t). We see that the dynamics of the density matrix p{t) drops out of the final formula, 
as does the interaction potential in the Hamiltonian H, leaving a result that is just the sum of 
contributions from the kinetic terms of the individual particles. 

We give now several specific examples of the formula (]7ip . First of all, in the standard white 
noise CSL model, one has uncoupled space and time correlators of the product form 

D{x - y,t - s) = G{x - y)6{t - s) , (72) 

which taking account of the fact that Jq ds6{t — s) = 1/2, gives 

F{x-y,t) = ^G{x-y) . (73) 

The spatial correlation function G{x — y) is the autoconvolution of the function g{x) introduced 
as the CSL smearing function. 



G{x - y) = I d zg{x - z)g{y - z) , 

) 



g{x)=(V2^rc\ e-^"'/(2r2,) , 



(74) 



which, incidentally, gives an alternative form of factor decomposition for this model. We will 
continue, however, to use the factor decomposition given by the Fourier transform, which is 



(2^ 



G{S-y)= I j—^cos[k.{x-y))e-^'^c , (75) 



so that substituting (j75p and (j73p into (j70p and writing k = \k\ we have 



F{k,t) = \e~''^'^ . (76) 



Substituting this into ()7ip gives for the white noise CSL model 



87 ^-^ mf 3A ^-^ m- 



2-^ M, A'm'^.rl. 2^ 



' 2,2-K^/'^rl; ^^ Mi Am%rl^-^M^ ' 



(77) 
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with TTijv the nucleon mass and A = 7m^/(87r^'^r^) the alternative form of the noise couphng 
generaUy used in the CSL hterature.^ This result agrees with the standard answer for the constant 
energy production rate in the CSL model. 

Consider next a variant of the product correlator model, in which there is a cutoff in the 
frequency spectrum, obtained by replacing j5{t — s) in (I72p by 

^'y{uj)it — s) = — / duj'y{ij) COS iuj{t — s)] . (78) 

This replacement turns the original coupling 7 into a frequency dependent coupling 7(w), with the 
specialization back to constant 7 given by (5^(tj)=-y = 7(5(t — s). In this case we find that 

/ ds5^{uj){t - s) = - 7(-) — sinu , (79) 

JO TT Jo ^ f^ ^ ^ 

which approaches the constant 7(0)/2 as t ^ 00, with the entire contribution in the infinite time 
limit coming from the infrared region of the integral near a; = 0. Thus even with a high frequency 
cutoff, there is a constant energy production rate at large times in a model with uncoupled space 
and time correlators. To avoid getting a constant energy production rate in the product correlator 
model, one must include an infrared cutoff, by taking 7(0) = 0. 

Finally, anticipating our discussion below of thermal noise, consider a correlator of the general 
form 

D{x-y,t-s)= ——^ — fik)cos(%- {x-y)j cos [uJk{t- s)j , (80) 

with LOk a wave number dependent angular frequency. Integrating to form F{x — y,t), we have 

F{x-y,t)= dsD{x - y,t - s) = / ^ f{k) cos (k ■ {x - y)) sm{ujkt) , (81) 

Jo J (2vr)^w^ V y 

which identifies the Fourier transform as 

F{k, t) = A/ sin(a;fct) . (82) 

Substituting this into ([7T]) gives for the energy production rate 

d^ , , sr^ mf f d'^k k'^f(k) , , , , 



^ In the CSL model literature, what we here call 'Ym% is termed 7, because the noise there is introduced as coupled to 
the nucleon number density rather than the mass density. Also, we note that the dimensionality of 7 is determined 
by the dimensionality assigned to the field (p, and is not the same in our white noise and thermal model examples. 
In the white noise CSL model, what we call 7 has dimensionality mass~^ in microscopic units with h = c — 1, 
whereas in the thermal noise model discussed below, where <j!> is taken as a conventional dimension one boson field, 
7 has dimensionality mass^^ in microscopic units. 
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Even when w^ oc A;, this expression is strongly convergent in the infrared as a consequence of the 
vanishing of phase space for small k values. Hence if f{k) is cut off sharply at large k values, as 
expected in thermal models, it leads to a vanishing energy production rate at large times by use 
of the Riemann-Lebesgue theorem. Integrating to find the total energy production ATrHp{t) = 
Ti Hp{t) - TtHp{0), we find 

H JO ^k 

which as t ^ oo gives, again by an application of the Riemann-Lebesgue theorem. 



ATVifp(t) = ^J]^jJ|"dfc^^[l-cos(u;fet)] , (84) 



A-K»)^iEif -^ ^ (-) 



B. Gamma radiation from atoms 

An important constraint on noise model parameters is provided by the spontaneous emission of 
gamma rays from atoms, a process first calculated for free electrons by Fu [6] and later calculated 
for general atomic systems by Adler and Ramazanoglu [71]. Results were given in the latter paper 
for a correlator of the form G{x — y')^-^!^) (t — s). Remembering that the CSL definition of 7 is mj^ 
times the definition of 7 used in this paper, and comparing ([73|) and ([75]) with ([80|) . we see that 



the results of 7| can be converted to apply to a correlator of the form of (fSTI) by the substitution 

7(w) _I,2„2 f(k) ^, , .. 

n 

When uJk has the form Wfc = \/fc2 + p?, making this substitution into (44) of [7] gives as the formula 
for the power radiation dP per unit photon energy dp from a hydrogen atom. 



dp 



1 
1 



[1 + (pao/2)2]2 






(87) 



with e^ ~ 1/137.04 and with ao = l/(e mdcctron) = 0.529 x 10 ^cm the Bohr radius. 

VI. MODELS FOR THE CORRELATION FUNCTION 

We turn in this section to a discussion of specific models for the correlation function D{x — y, ti — 
t2) introduced in ([5]). We first briefly consider the standard CSL factorizable correlation function 
with white noise, and its variant with a cutoff in the noise spectrum, which has been the basis of 
most discussions to date of the phenomenology of objective reduction models. However, one would 
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in general expect the spatial and temporal structures of the correlation function to be intertwined, 
and in particular, a correlation function arising from fields with a particle interpretation will have 
the spatial and temporal correlations coupled by a mass-shell constraint. This is the motivation 
for the models discussed in the remainder of this section, which are based on a classical model 
extracted from the quantum thermal Green's function of a boson of mass /i. 

A. The product correlator model 

The product model for the correlation function was written down above in ()72|) through ()79p . 
With a white noise spectrum, the standard CSL choice for the noise strength parameter is m'j^'y = 
lO^'^'^cm'^s"^, and the standard choice for the correlation length is re ~ lO^^cm. The white noise 
model with these parameter choices obeys all experimental upper bound constraints, and readily 
explains measurements in which riout ~ 10^^ nucleons are displaced by a distance of at least re- 
in [a| Adler gave a reanalysis of the upper and lower bounds on parameters in stochastic reduc- 
tion models. Under the assumption that latent image formation, in either photography or etched 
track detectors, constitutes a measurement (rather than the measurement occurring only through 
the subsequent development that reveals the latent image), he concluded that the noise strength 
parameter 7 should be larger than conventionally assumed in the CSL white noise model, by a 
: 'actor of 2 x 10^ . This however conflicts with bounds set by Fu [6] and Adler and Ramazanoglu 



a 



on spontaneous 11 keV gamma radiation emission from germanium, unless the white noise 
spectrum is cut off at energies below 11 keV by the spectral weight 7(0;) appearing in ([78]) . Such 
a cutoff would still allow sufficiently rapid state vector reduction to account for observed measure- 
ment times, as already noted in the review of Bassi and Ghirardi [9]. Thus, the product model for 
the correlation function, with a high frequency cutoff in the noise spectrum, is consistent both with 
all upper bounds, and with the assumption that latent image formation constitutes a measurement 
signaling state vector reduction. Such a correlation function rnight arise from a pre-quantum theory 



in which quantized fields are not the primary entities, as in 



10|. But as already noted, a product 



correlation function is not expected to arise from quantum fields with a particle interpretation. 

B. Thermal correlation function model 

In this section we shall motivate a model for the correlation function D(x — y,ti — ^2) by 
considering the correlation function for a quantum field in a thermal state at temperature T. Let 
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(x, t) be a scalar quantum field, with the mode decomposition 



{x,t) 



d^k 



1 



1/2 



a(A?)e*(^-^-^^-*) + at(A?)e~*(^-^~^'=*) 



(88) 



-2wfe(27r)3. 

where a(A;) and a^(A:) are the mode annihilation and creation operators, and where the mode energy 
Wfc is 



i^k 



fc2 + ^2 ^ 



(89) 



with fi the scalar field mass. We have, as before, set Planck's constant h equal to unity, and 
also set the Boltzmann constant equal to unity, so that in a thermal state at temperature T the 
expectations of products of creation and annihilation operators are given by 

{a{k)aHk'))=5^{k-k')[l + N{k)] , 
{a^k)a{k'))=6^{k-k')N{k) , 



(90) 



with the mean occupation number N(k) given by 



N{k) 



CT — I 

From these equations we can now calculate the correlation function 



(91) 



(0(f,ti)0(y,t2))= / d-'k 



2a;fc(27r) 



[1 + Nik)y 



+ N{k)e~ 



(92) 



which can be written as the sum of a temperature-independent part A-^-{x — y,ti — t2) and a 
temperature-dependent part D{x — y,ti — ^2), as follows, 

{(j){x, ti)4>{y, t2)) =A+(x -y,ti- t2) + D{x - y, ti - t2) , 



D{x-y,ti -12) 



Sk 



1 



2wfe(27r)2 



Sk- ^(^) 



2wfe(27r) 



iyk-{x-y)-uik{ti-t2)) _|_ ^-i[k-{x-y)-ujk{ti-t2)) 



(93) 



In the zero temperature limit, D{x,t) vanishes, and (j93p reduces to the temperature-independent 
piece A_|_, which is one of the standard relativistic quantum theory vacuum Green's functions 
arising directly from the non-commutativity of a{k) and a'{k), and is a complex number for general 
arguments. The real-valued temperature-dependent piece D{x — y,ti — ^2), on the other hand, is 
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invariant under the interchange x, ti ^^ y,t2, and therefore can serve as a model for the expectation 
of real, classical, commuting noise fields introduced in ([5]). 

Since A^(A:) and cj/c are even in A;, writing e *^ '^^~y)) = cos (fc • (x — y)) i^sin (fc • (x — y )), the 
sine functions average to zero, and the formula for D(x — y,ti — ^2) simplifies to 

/d'^k N(k) 
■7—3 cos(A;- (f -y))cos(u;A:(ti - ^2)) , (94) 

which has the form assumed in (j80p . with f{k) = N{k) as given in (j9ip and with uJk the energy- 
momentum relation given in ()89p . We shall slightly generalize the model specified by (I94p and ()9ip . 
by introducing a thermodynamic chemical potential C into the occupation number, which we thus 
write as 

N(k) = ^ , (95) 

e T — 1 



which allows us to accommodate systems with general particle density [111 ]. For the case of noise 
fields associated with particles having a standard energy- momentum dispersion relation, ()94p . (I89p . 
and (f95]l constitute our basic model for the correlation function D{x — y,ti — 12). Corresponding 
to this model, the function F{x — y, t) defined in (jlip is given by 

/* f 6?k Nik) 

dsD{x - y,t - s) = -—^ — —cos{k-{x-y))sm{uJkt) , (96) 

and the integral appearing in the rate function r(t) of (J17p is given by 

I{x-y,t)= dsF{x-y,s)= I ——^ — 3- cos (fe • (f - y )) [l - cos (w^t)] . (97) 
Jo J i^^j ^k 

C. Dilute and nonrelativistic limits 

Let us consider now the dilute limit of (j94p and (j95p . obtained [ill] by letting the chemical 
potential C, be large and negative, so that N{k) becomes 

Nik) ~ e~^ . (98) 

We will be particularly interested in applying (j98]) to the nonrelativistic case T « /j,, where we 

can expand 

fZ p 

UJk = y k'^ + fi^ c^ fi + — , (99) 

so that N{k) becomes 

iV(fe) = e-('^-^)/^e-^''/(2MT) (lOQ^ 
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Where w^ appears as a denominator factor in (|94j) . (j96|) . and (j97j) . it can be approximated by //, 
so these equations become respectively 



D{x,t) 



F{x,t) 



I{x,t) 



^I'-O/T r d^k 



M 



H- 



(27r)3 



'if'-O/T f ^3^ 



(27r)3 



/^^ 



-(/^-^/^ /■ d^k 



(27r)3 



^fe2/(2MT) 



-fc2/(2MT) 



-kV(2f,T) 



k^ 



cos (A; • x) cos ((^ H )t) , 

^ ^ 2;U 

cos (A; • x) sin ((/i + —)t) , 
cos (A; • x) [l — cos ((/u + — )i)] 



(101) 



Carrying out the angular averaging over k, remembering that it is the difference -D(0, t) — D(x,t) 
that enters into the reduction formalism, and writing R= \x\, k = \k\, (jlOip yields 



D{6,t) - D{x,t) 



F{6,t) - F{x,t) 



/(0,t)-/(f,t) 



(^-f)/^ /-"^ k^dk 



M 



e-(^-^)/^ z-'- edk_ 

' f^' y-oo(2vr)2' 



^l- 



-ky{2^,T) 



-kV(2tiT) 



,-ky(2^^T) 



sm{kR) 



kR 
sm{kR) 



kR 
sin{kR) 



kR 



^2 

COS {in + —)t) , 

k^ 
sin((/i + — )i) , 

[l-cos((^ + ^)t)] . 



(102) 



The integrals in (jl02p can all be evaluated from the formula 



2 , 2^sin(a;/3) 

x dxexp[—ax )- 



vr 



-/3V(4a) 



(103) 



x(3 2a3/2 

with results that are summarized in Appendix B. In particular, for large times, the formula for 
/(O, t) - I{x, t) limits to 



10, t = oo) - I{x,t = oo) ^ e-^^-'^y^(—Y^' \l - e-«'^^/2 

\2TTflJ 



(104) 



and the formulas of Appendix B show that the characteristic reduction time tji for approach to 
the asymptotic value of (jl04p is the inverse temperature T~^. 

To compare this to the standard CSL model formulas, let us look at the decay of the off-diagonal 
density matrix element of a one particle system of mass equal to the nucleon mass rriN, which we 
have seen is governed by the same rate function T as state vector reduction. From ()16p and (J19p . 
at large times one has 



{f^\p{t = oc)\f^) = e-r(°°)(f V(0)|f 



^2\ 



(105) 
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with 



r(t = oo) =2jmj^[I{0, oo) - /(f ^ - f ^ oo)] 

VzTr/i/ 



(106) 



where we have written R = |r ^ — r ^|. The comparable formula in the CSL model is given in (8.15) 
ofU, 

(fXt)|f2) = e-r^''W(fi|p(0)|f2) , (107) 

where r*-"^'^(t) is given by 

and where (as remarked above in a footnote) j*^^^ is what we call ^mj^. We see that the functional 
form of the -R-dependence in (jl04|) and (jlOSp is the same, with the CSL model correlation length 
re related to the nonrelativistic thermal model parameters by 

9 1 / T \3/2 / 1 x3/2 

'■^=2]^ • (2;^) =" (S^) ^ <™' 

However, whereas F (t) grows linearly with time for large times t, in the thermal noise model 
r(t = 00) approaches a constant. This means that to achieve the degree of density matrix diago- 
nalization, or state vector reduction, attained in the CSL model in time At, the parameters in the 
thermal model must obey 



^^^CSL ^ ^7^g-{;.-C)/T (110) 



D. Can thermalized dark matter be the noise source? 

As we have already noted, one motivation for studying non- white noise is to investigate whether 
there can be a cosmological origin for the noise that drives state vector reduction in objective 
reduction models. Since there is now strong evidence that about a quarter of the closure density 
of the universe consists of dark matter, and since weakly interacting massive particle (WIMP) 
candidates for dark matter are expected to be thermalized, it is natural to apply the results of the 
preceding section to an analysis of whether dark matter can account for the noise coupling in ([6]) . 
We will not attempt to discuss here the necessary conditions for dark matter to give a real-valued, 
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as opposed to an imaginary- valued, noise term in the Schrodinger equation; this important question 
will be deferred to future work. What we shall do in this section is to assume that a real-valued 
noise coupling can be achieved, and to investigate the phenomenological implications of assuming 
that state vector reduction is associated with observed dark matter parameters. 

A few basic facts about dark matter are needed. If dark matter is due to WIMPs, then observa- 
tional evidence [12] suggests a WIMP distribution in the galactic halo of mass density Pmass = 0.3 
GeV cm~^, and a Maxwellian velocity distribution with Uj-ms = 220km s~^ = 7.3 x 10~^c. The 
r.m.s. velocity is estimated from the formula 

which describes the gravitational binding of WIMPs of mass // to the galaxy of mass Mgaiaxy, at 
radius rgaiaxyi with G the Newton gravitational constant. Direct limits on possible solar system- 
bound dark matter are weaker [13l |. [14| by a factor of 3 x 10^, that is, pmass ss ^ 0.9 x 10^ GeV 
cm^'^. There is at present no observational limit on possible earth-bound dark matter. If there 
were solar system-bound dark matter, around the radius of the earth's orbit the r.m.s. velocity, by 
(jllip would be Urms ~ 30km s^^ = 10~^c, and for earth-bound dark matter, at the radius of the 
earth's surface, the r.m.s. velocity would be frms ~ 8km s~^ = 0.27 x 10~^c. 

Because the WIMP mass /i cancels out of (11111) . there is currently no direct information about 
the dark matter particle mass. Dark matter particles coupling to the mass density cannot be 
too light, or they would conflict with gravitational fifth force experiments. If we write the noise 
coupling as 

'' = w • ("^' 

then the fifth force experiments require 

exp(-M/^5) ^ 1 .„„,, 

IVl ^^-'Planck 

with ;U5 the fifth force scale limit, currently [15|] around ^5 ~ 1.4 x 10^'^ eV. This gives the lower 
bound on M, 

^j> ;LQl9-0.22M/M5QeV . (114) 

In addition to this constraint, there are also mo del- dependent astrophysical limits on the dark 
matter mass; for example, warm dark matter candidates must have masses greater than 1 keV |16l |. 
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For a Maxwellian distribution with N{k) given by (jlOOp . the r.m.s. velocity is given by 



2 3r 



^rms = , (115) 



SO that using (jlOOp we have 



v3/2 

frms = • (116) 

Ijrc 



Hence for a given dark matter r.m.s. velocity, the correlation length re and the dark matter 
temperature T are determined as functions of the dark matter mass /i, 

/372 



re 
T 



r'^rms 



(117) 



Integrating N{k) over phase space, the number density pn is given by 

/d^k 
7^:33^(^) 



(27r)3 87r3/2r3, 



which determines the factor containing the chemical potential C in terms of pm, p and re, 



(118) 



g-(M-C)/T ^ ^8^3/2^^ . (119) 

From these equations, together with (jl06p and (jl09p . and the assumption that the lower bound 
of ([55]) gives a good approximation to the reduction factor,^ we get the following estimate. 



Reduction factor ~e-2'^(*=°°) , 

(120) 



2r(t = oo)=4(— j -^niV 



^ Note, however, that although r(t = oo) is positive, the uniform positivity assumption on the integrand used to 
derive the upper and lower bounds is not obeyed in the thermal model; see the formulas in Appendix B. Also, the 
simple model analyzed in Appendix D gives exponential reduction as in the lower bound, but with a reduction 
factor e~'"'*^°°\ and so the rates calculated from the lower bound may be optimistic by a factor of two. 
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Here, in accordance with the properties of T discussed in Sec. 2, n is the number of displaced 
nucleons that are bunched within a correlation length re, and N is the number of such bunches of 
displaced nucleons. 

Using (J117p and (jl20p . we can now make some estimates of the effectiveness of thermal dark 
matter in producing state vector reduction in the mass-density coupled model. Rewriting ()117p in 
the form 

0.24 X 10-13 iQgV 



re 



--tR 



cm 

_ 0.2 X 10-23 iGeV^ 

~ 2 ' 



(121) 



we get the following tables of values. For the correlation length re in the body of the table in 
cm, versus the dark matter mass /u in keV and its r.m.s. velocity appropriate to the galactic halo 
{vh = 220kms-i), solar system-bound dark matter {vs = SOkms-^), and earth-bound dark matter 
{ve = Skms-"*^), we have 



M^ 


1 


10 


102 


103 


10" 


10^ 


Vh 


3 X 10-5 


3 X IQ-e 


3 X 10-^ 


3 X 10-^ 


3 X 10-9 


3 X IQ-" 


Vs 


2 X 10-4 


2 X IQ-s 


2 X 10-s 


2 X 10-^ 


2 X 10-8 


2 X IQ-i" 


Ve 


9 X 10-'' 


9 X 10-5 


9 X 10-6 


9 X 10-^ 


9 X 10-8 


9 X IQ-i" 



TABLE I: Correlation length re (cm) versus /j, (keV) and r.m.s. velocity v. 

Similarly, for the reduction time tR in seconds in the body of the table, versus the dark matter 
mass and its r.m.s. velocity, we have 



M-^ 


1 


10 


102 


103 


10" 




10*^ 


Vh 


4x 10-12 


4x 10-" 


4 X 10-" 


4 X 10-15 


4x 10- 


-16 


4 X 10-18 


Vs 


2 X 10-1" 


2 X 10-" 


2 X 10-12 


2 X 10-13 


2 X 10" 


-14 


2 X 10-16 


Ve 


3 X 10-9 


3 X 10-1" 


3 X 10-11 


3x10-12 


3x 10- 


-13 


3 X 10-15 



TABLE II: Reduction time tn (s) versus /i (keV) and r.m.s. velocity v. 



Solving (jl20p for the value of jpm which yields 2r(t = oo) = 1, which is the minimum value of 
the exponent beyond which reduction of the state vector starts to occur, we get 

1.5 X 10^3 / ^ 
l^lGeV 



IPn 



n 



2N 



GeVcm 



-1 



(122) 



27 



From this, we get further tables of values. For ^Pm in the body of the table, in GeVcm"^, versus 
the dark matter mass p. in keV, and the effective number of displaced nucleons riout = n^N = 10^^ 
corresponding [9| to the standard CSL model, or riout = n^N = 10^ corresponding to estimates 
based^ on latent image formation, we have 



M^ 


1 


10 


102 


10^ 


10^ 


10^ 


1022 
108 


2 X 10-21 
2 X 10-^ 


2 X 10-13 
2 X 10-5 


2 X 10-1^ 
2 X 10-3 


2 X 10-15 
2 X 10-1 


2 X 10-13 
2 X 10 


2 X 10-9 
2 X 10^ 



TABLE III: ^p^ (GeV cm i) versus p (keV) and Uont- 

If we make the assumption that 7 = l(TeV)-^ = lO-^(GeV)-^, we get a table of values giving 
Pm in the body of the table, in GeVcm-^, versus the dark matter mass and the effective number 
riout of displaced nucleons. 



p^ 


1 


10 


102 


10^ 


104 


10*^ 


1022 
108 


3 

3 X lOi'' 


3x 10* 
3 X lOi** 


3x 10^ 
3 X 1022 


3 X IO12 
3 X 1026 


3 X IOI6 
3 X 10^0 


3 X 1024 
3 X 10^8 



TABLE IV: pm (GeV cm ^) versus p (keV) and nout- 



From these tables, we see that state vector reduction, by the standard CSL criterion (riout = 
10^^), and with a correlation length within a decade of the standard CSL value re = IQ-^cm, 
is achievable in the dark matter model for dark matter masses in the range of 1 to 10 kilovolts, 
with 7 ~ ITeV- and with pm below the current upper limit on solar system-bound dark matter. 
Adopting the latent image criterion (riout = 10^) requires dark matter densities that are much too 
large, so either the latent image analysis of [8j] needs modification, or the dark matter model is 
unworkable. 

For a dark matter mass ^ of a kilovolt or greater, and the current limit on the fifth force scale 
//5, the fifth force bound of (J113p becomes 



M> 10 



19-0. 15x10*5 



(123) 



■1 In the CSL model, one assumes n = 10®, which is the number of nucleons in a volume of linear dimension 10 ^cm, 
and iV = 10", giving n?N = 10^^. The latent image estimates of [1] take n = 5640 and iV = 20, giving n^N ~ 10*. 
The CSL model assumes a reduction rate of lO^s-^, whereas the latent image estimates assume a much smaller 
reduction rate of 30 s~^ , which is why in a white noise model the ratio of the noise strengths between the two 
cases is ~ 10®, rather than the ratio ~ 10^* of the ri^N values. 
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which is strongly obeyed for the M values in the GeV to TeV range that are interesting. Referring 



to the discussion following ([94]) . and using ([95]) and ()117p . we see that the function /(\/p^ — fjfl) 
in the formula (j87p for the radiated gamma power from a hydrogen atom becomes 

f^/„2 _ ..2\ - ^ p-(M-C)/T -(p-A*)/?^ 



(124) 



Since for /i in the 1 to 10 kilovolt range and for p = 11 kilovolts, we have 



T,Hp{t = oo) = ''•'^^ ^7 . (126) 



^^^^>6xl05 , (125) 

the negative exponential in the final factor of (jl24p dominates all other factors in this equation 
and in (I87p . and so the experimental bound on 11 keV gamma radiation is strongly satisfied. 

For both values of nout displayed in the tables, the reduction time is sufficiently rapid, shorter 
than a few times 10~^ seconds, to account for realizable measurements. Finally, the total energy 
imparted by the noise to an isolated nucleon is obtained by evaluating ()85p by using the form for 
f{k) in the dilute nonrelativistic thermal model, giving 

For the CSL value of nout, this is smaller than 10^^^ degrees Kelvin for all values of the dark matter 
velocity and mass in the tables, and so is sufficiently small so as to be unobservable. 

The conclusion from this analysis is that, i/dark matter couplings to ordinary matter have the 
anti-self-adjoint component needed to give a real- valued noise term in the Schrodinger equation, and 
if dark matter densities in the vicinity of earth are larger than the galactic halo density, but within 
current limits on solar system-bound dark matter, one could realize the standard CSL reduction 
model with the standard parameter values, and obey various important experimental constraints. 
The italicized assumptions make this mechanism for realizing state vector reduction conjectural; at 
worst, we have given an interesting toy model for reduction incorporating a non- white noise with 
a mass-shell constraint. 



E. Thermal unparticles as the noise source 



Recently Georgi 



231] has introduced the concept of what he terms an "unparticle" , a field 



characterizing a scale-invariant sector of a low-energy effective field theory. This is of interest for 
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collapse models, since if the noise field of Q is the low-energy manifestation of a pre-quantum 
dynamics, such as discussed in the book pi|, it is plausible that it could have a scale-invariant 
structure. Moreover, such an unparticle field, if a cosmological relic field, will have a thermal 
correlation structure. The concept of thermal unparticles has been introduced in a recent paper of 
Chen et al. , 2j], who construct the thermal unparticle partition function by using the observation of 
that an unparticle field can be constructed as a field with a continuous distribution 



3, 



Krasnikov 

of mass /i^, characterized by a scale invariant spectral function p(/i^) oc {ij?)'^~'^ . More specifically, 
one obtains the unparticle propagator and partition function by integrating the corresponding 
propagator and partition function for a scalar field of squared mass ^"^ over the range < /i^ < cxd, 
with weighting function p{^J?') = {d— l)h?^^^'^' {^'^Y^'^ . Here d is the anomalous scaling dimension 
characterizing unparticle physics, and A is a scale parameter (the cutoff for the low-energy effective 
theory) with dimension of mass.^ 

In Appendix E we use the same method to construct the unparticle thermal correlation function 
from the thermal correlation function of (j94p and (j95p for a scalar field of mass ^■^. From this 
correlation function, we calculate the integrals needed to study both the state vector reduction 
rate and the noise-induced energy production. We recapitulate here two key formulas obtained 
from Appendix E, both of which apply to a one particle system of mass m. For the decay rate 
r(t) of the off-diagonal matrix element (x|/7(t)|0), which we have seen is also the reduction rate, 
we have 

T{t) =2W[Iu{^,t) - Iu{x,t)] 

=- 2 / ^ ^^^c —^ / dv\V-co-s{v^\x'\)\(\-v^f-^ , 

^ h e~^r- — 1 Jo 

(127) 

where the subscript U on I corresponds to the notation of (jE16p of Appendix E. For the noise- 
induced energy acquisition rate and total energy acquired by a particle of mass m, we have from 
([E20l) and (|E2T]) of Appendix E, 

* (2ir)2 r(3/2 + d} J„ eT^ - 1 



^ Strictly speaking, the integration over fi^ should extend only up to A'^, but when the temperature T << A, the 
integration for the partition function and thermal correlation function is effectively cut off by N{k) of (|95p . and 
so negligible error is made in extending the upper limit to oo. 
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and 

Turning our attention first to (jl27p , we note that the inner integral over v is always convergent 
at f = 0, and is convergent at i; = 1 for Red > 0. Because the inner integral in p27|) vanishes 
as w^ for small uj, the integral over uj in (jl27p has precisely the same convergence properties at 
uj = as the integral giving the total energy production in (|129p . To study convergence, there are 
two cases to consider, (i) the chemical potential C is negative and nonzero, and (ii) the chemical 
potential ( is zero.^ 

In the first case, of strictly negative C, the denominator e t — 1 is nonzero even at a; = 0, and 
the integrals of (jl27p and (jl29p converge at a; = even when the factor 1 — cos (wt) is replaced by 
unity, as long as Red > 0. So in this case we can extract the infinite time limit by invoking the 
Riemann-Lebesgue theorem, and simply dropping the term cos (tot) in (|127p and (jl29p . giving the 
formulas 

/ du 

e""T" -I Jo 

and 



r(oo) = ^^ 2 / dio-:;^^ / dv[l-cos{vuj\x\)]{l-v^)'^~'^ , (130) 

^ Jo p.^r- — 1 Jo 



(2i)2 r(3/2 + d) y„ e=i* _ 1 ^ ' 

Corresponding to the fact that the total energy production is finite, the energy production rate of 
(|128p vanishes at large time. Referring now to (|130p . we see that there are two subcases governing 
the large |x| behavior, which we call (ia) and (ib). In subcase (ia), corresponding to Red > 1, 
the CO integral is convergent without using the u;^ factor arising from the inner integral. So in this 
subcase we can apply the Riemann-Lebesgue theorem to the inner integral in the limit of large \x\, 
by dropping the term cos (wu;|x|), leading to the conclusion that r(oo) varies from at |x| = to 
a finite value at |x| = oo. In subcase (ib), corresponding to 1 > Red > 0, the o;^ factor from the 
inner integral is needed for convergence, and on changing integration variable from uj to u = uj\x\ 
one sees that r(oo) grows as |xp'^~'^) as \x\ -^ oo. 

In the second case, of vanishing chemical potential C, the denominator e t — 1 vanishes at 
a; = 0, and the integrals of (jl27p and (jl29p now behave for small cj as 

^j^2rp^2{l-d) I- j-l 

T{t) ~ 2 / duu'^^-^ [1 - COS (cjt)] / dv[l- cos (vcj|x|)] (1 - v^f'^ , (132) 



The chemical potential must always be less than or equal to zero, so there is not a third case of positive (", which 
would correspond to a physical region pole in the integrands coming from the vanishing of the denominator e ^ — 1 
in all three integrals. 
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and 

T.^„*) - T.Hm ^ '-"^'^"^^^ I ^--' [1 — (..)] . (133) 

There are now two subcases, which we label (iia) and (iib). In subcase (iia) we have d > 1/2, and 
both integrals (1132p and (jlSSp converge at w = without using the uj'^ factor that comes from 
1 — cos (wt) . So in this case, which behaves much like case (i), we can apply the Riemann-Lebesgue 
theorem to take the limit as t — > cxd by dropping the term cos (wt) , leading to finite values for r(oo) 
and TrHp{oo) — TiHp{0). One can then proceed to analyze the large \x\ behavior of r(oo), as was 
done previously in case (i), with the conclusion that this is finite for d > 3/2 and it behaves as 
|^|3-2d £qj, 2^2 > d > 0. In subcase (iib), we have 1/2 > d > 0, and the uj"^ coming from the factor 
1 — cos (tijt) is needed for converence; defining a new integration variable u = cot, we see that both 
r(t) and TrHp{t) grow as t^^'^'^ in the large t limit, and correspondingly, the energy production rate 
decreases as t~'^'^. So for vanishing chemical potential, and 1/2 > d > 0, we have the interesting 
situation that one achieves perfect reduction at infinite time (that is, r(cxD) = cxd), although the 
reduction rate and the total energy production both grow as a fractional power of t, rather than 
linearly with t as in the standard CSL model. Correspondingly, the energy production rate vanishes 
as a fractional power of t at large time, which should make it easy to satisfy cosmological constraints 
81] on the noise strength parameter. 

We conclude that the thermal unparticle model exhibits a range of interesting behaviors, de- 
pending on the values of the chemical potential C and the unparticle dimension d. In addition to 
these two parameters, the effective noise strength ^A^*^^~^ and the temperature T are also param- 
eters of the model. Given the complexities of this four-dimensional parameter space, we do not 
attempt phenomenological fits of the model to experimental constraints on the noise strength, but 
this is clearly an interesting topic for future investigation. 

VII. SUMMARY AND DISCUSSION 

We now summarize what has been done in the preceding sections and what is in the appendices, 
and sketch some directions for extensions of our investigations. In Sees. I-V we have continued 
the study of non-white noise models initiated in (I), focusing on the special case in which the noise 
field couples to the particle density. The analyses of Sees. II, III, and Appendix D identify the 
characteristic rate functions governing density matrix diagonalization and state vector reduction, 
and show that both processes are exponential with the same rate function, in the simplified case 
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(a single particle in a superposition of two localized states) discussed in Appendix D. In Sec. IV, 
we completed our formal analysis for non-white noise by deriving the corresponding Fokker-Planck 
equation, allowing us to make contact with earlier work of Pearle [5|]. In Sec. V, with an eye 
towards phenomenological applications, we analyzed energy production and gamma radiation by 
atoms in terms of the correlation functions of the non-white noise model. 

In Sec. VI we turned to a discussion of specific models for the noise correlation function. After 
a brief discussion of the product correlator model that has been the basis of most earlier work on 
objective state vector reduction, we turned to a detailed analysis of a thermal correlation function 
model, in which the spatial and temporal correlations are linked by a mass-shell constraint. We 
showed that the dilute, nonrelativistic limit of the thermal correlator model can be put in direct 
correspondence with the formulas of the standard Gaussian CSL model. We then gave a detailed 
phenomenological analysis of thermal dark matter as the noise source, and concluded Sec. VI 
with a discussion of the behavior of thermal unparticles as the noise source, sketching qualitative 
behaviors for a range of values of the chemical potential and of the unparticle anomalous scaling 
dimension. The examples given included cases in which T{t) and the energy production both are 
finite at t = oo, and in which T{t) and the energy production both grow as a fractional power 
smaller than unity as t ^ oo. 

The appendices deal with various details connected with the main discussion. In Appendix A 
we estimate the validity of the Markovian approximation used in the energy production discussion, 
while in Appendix B we compare the master equation used in our discussion with a more general 
class of master equations appearing in the literature. Appendix C gives the evaluation of integrals 
for the dilute, nonrelativistic model, while Appendix E gives details of the unparticle correlation 
functions. Appendix D shows that, in a simple model, reduction is exponential in the rate function 
r(t), indicating that the lower bound derived in Sec. Ill, as opposed to the upper bound derived 
there, gives the better estimate of the qualitative reduction behavior. 

We can point to a number of possible directions for generalization or extension of the results 
of this paper, (i) We have considered only the case of a real noise coupling, corresponding to an 
anti-self-adjoint Hamiltonian term. More generally, one could consider a complex noise coupling, 
containing both real and imaginary noise couplings, with the real term contributing both to density 
matrix diagonalization (i.e., decoherence) and to state vector reduction, and the imaginary part 
contributing only to decoherence. (ii) For simplicity, we have only considered a scalar noise field 
(j), but a general treatment of non-white noises would allow for the possibility of spin-1/2 or spin- 
1 noise fields. Such an extension may ultimately be required on phenomenological grounds to 
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make contact with experiment, (iii) Although we sketched the quahtative behavior of the thermal 
unparticle case, we did not attempt to make a quantitative phenomenological survey of the four- 
dimensional parameter space of this model, and this would be of interest, (iv) The derivation of 
lower and upper bounds on the reduction rate in Sec. Ill made use of a positivity assumption, 
which is not obeyed in the thermal correlator model; can this assumption be eliminated? (v) The 
model calculation of Appendix D indicated an exponential dependence of the reduction factor on F, 
agreeing with the corresponding density matrix diagonalization calculation but differing by a factor 
of 2 in the exponent from the corresponding lower bound of Sec. III. Can this result be generalized 
to the case of many particles and a wave function that is the superposition of many localized states 
as in (j3ip ? Clearly, a general argument that the reduction factor has exponential rather than power 
law dependence on T would be significant for the phenomenology of objective reduction models. 
(vi) In Sec. VIB we formulated our thermal model for the correlation function, by neglecting the 
temperature-independent Greens function A^-, which reflects the non-commutativity of creation 
and annihilation operators. As noted, this gives an effectively classical model for the thermal noise. 
It would be worth exploring a fully quantum-mechanical treatment of state vector reduction by 
a thermal noise field, in which all parts of the quantum mechanical correlation function (j92p are 
retained. 
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APPENDIX A: MARKOVIAN APPROXIMATION 

One can estimate the validity of the Markovian approximation by considering the case of a 
single free particle of mass m, so that H = p^ /{2m). Then we easily calculate that 

M{y,s- t) = m5^{y-e'"^''*Ue"'"^'~*^) =m5^{y-q- {p/m){s - t)) , (Al) 

so that repeating the steps leading to ([83]) we find 

^TTpit)H = ^m fds [^^ cos Uit-s))T:tOis-t) , (A2) 

at Jo J (27r)-^ Wfc V / 

with 0{s — t) given by 

0{s -t) = _i[e-^^-(9+(p7™)(^-t)) ^ [e^^'-9,p2]] . (A3) 

This expression can be simplified by use of the Baker-Hausdorff theorem and the canonical com- 
mutation relations, giving after considerable algebra, and dropping terms that are odd in k, 

Ois-t)=ecos(^is-t))cos(^is-t))-2p.ksm[^{s-t))sm[^{s-t)) . (A4) 

We see that all dependence of 0(s — t) on s — t is through oscillatory terms. Assuming that the 
characteristic spatial variation scale of the problem is governed by Wfc ~ |A;| ~ |p | ~ /cmax, with fcmax 
the characteristic /c- value at which f{k) cuts off, then when the particle mass m is large enough 
for the kinetic energy at A^max to obey 

"^ << fcmax ) l^^j 

2m 

the variation oi 0{s — t) with s is much slower than that of the cosine factor in ()A2p . In this case 
the integral in (IA2[) is well approximated by replacing 0{s — t) by 0(0) = k'^, which recovers the 
result of the Markovian approximation made in Sec. V. 

APPENDIX B: COMPARISON WITH MASTER EQUATIONS FOR DECOHERENCE 

A further understanding of the effect of the thermal field (j){x, t) on the evolution of the wave 
function can be obtained by comparing ([5]) for the density matrix with typical master equations 
used for describing open quantum systems. Here we will follow the path outlined in [17|, where a 
coraparison of this kind has been made between the GRW model [18|] and collisional decoherence 



201, 



2li. 



191, 



35 



We consider the evolution of a single particle; under the Markovian approximation {M{y, s—t) 
M(y, 0) = M{y)) discussed in sec. V, ([8]) reads: 



d , , 



-i[H,p{t)]+Ct[p{t)], 



(Bl) 



with 



Cfip] 



J d^x d% [Mix),[Miy),p]] Fix-y,t), 



(B2) 



and M(x) = 'm6^{x — q). The term £[*, which includes the effect of the thermal field <j){x,t) on 
p(t), is the one we will focus on. Let us introduce the Fourier transform: 

d^k 



F{x-y,t) 



(2vr)^ 



F{k,t)e"'- 



ik-{x-y) 



(B3) 



with F{k,t) = F{—k,t) due to spatial inversion invariance. One can rewrite ()B2p in terms of 
F{k,t) as follows: 



Cfip] = 2m^j 



d^k 



F{k,t) 



^^k-qp^-ik-q 



(B4) 



(27r)3 

The above expression falls i nto the general class of translational-invariant Markovian master 
equations first given by Holevo 



22( 1 which, in the case of a bounded mapping £, reads: 



/OO p ^ 

dp{k)Y, e''^-^^Ln{k,p)pLiik,p)e-''^-'^ - - {^ (fc,p)L„(fe,p),p} 
n=l L 



(B5) 



where Ln{k,p) are bounded functions of the momentum operator p, and p{k) is a positive cr-finite 
measure. Briefly, the physical content of ()B5P is the following: the unitary operators e**"'"^ and e^'^^''^ 
induce a momentum transfer to the particle by an amount equal to k, while the operators Ln{k,p) 
imply that the momentum transfer to the particle depends on the momentum of the particle itself. 
This allows for mechanisms such as relaxation to take place. 

Equation ()B5P reduces to ()B4p under the following circumstances. Let us assume that Ln{k,p) = 
Ln{k) does not depend on the momentum p of the particle. They then become c- number functions, 
commuting with all other operators. By setting: 



d//(fc)^|L„(fc)|2 = 2m' 



7 



n=l 



d^k 

(2^)3 



F{k,t), 



(B6) 



the link is established. Of course, in the truly Markovian case one has D{x — y,t — s) = 
G{x — y)6{t — s) so that F{k,t) = (l/2)G{k) is independent of time, where G{k) is the 
Fourier transform of G{x — y). 
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According to the above analysis, the effect of the thermal field is not only that of localizing the 
wave function in space (this is a consequence of the specific form of the stochastic equation (I13p ). 
but also of exchanging momentum between the particle and the field. This is the reason why both 
the momentum and the energy of the particle are not conserved, in general. One would expect the 
energy of the particle to thermalize to that of the random field; however, the model described by 
(|13p does not allow for thermalization, since the operators Ln{k) do not depend on the momentum 
p of the particle. This is in agreement with the results of Sec. 5A on energy production. The 
comparison with decoherence suggests how the model can be modified in order to include also such 
an effect; this will be a subject of future research. 

APPENDIX C: INTEGRALS IN THE DILUTE, NONRELATIVISTIC THERMAL 

MODEL 



From (jl03p we find 



°° edk 



g ky{2^,T)^^}0^^ (f^y/^^-it.TR^'m 



oo 



(27r)2 kR V 27r / 



°° k^dk _fc2/(2^r)Sin(fci?) f,, , k^ 



exp u(/i + -— )t 



_o,(2vr)2- kR ^^^^V'- ' 2^ 

Ivr / 

exp (i{fit + (3/2) tan~^(tr) - {ptT'^R'^ /2)/{l + t'^T^] 



(CI) 
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from which, by forming hnear combinations, taking real and imaginary parts, and taking hmits as 

i? = |x| — > 0, we get 

D(6,t) - D(x,t) ^^ (^\^''^{l+fT'^Y^/^\cos (fit+{5/2)tan-\tT) 

jj \2tt J L \ 



COS i^t + (3/2) tan-^(ir) - {fitT'^R^/2)/{l + t^T^ 



F(0, t) - Fix, t) c^- 2 (^) (1 + i'T^)-'/' [ sin [^it + (3/2) tan'^ (tT) 

_g-(^TRV2)/(i+t^T2) g.^ / ^ ^ (3/2) tan-i(tr) - (/itr2i?2/2)/(i + ^2^2^ 



;i+t2^2)-3/4 



/(d,t)-/(f,0:.^:^(f)^^^{l-e-(^-V2) 



cosUt+(3/2)tan"^(tr) 
_g-(^TRV2)/(i+t^T2) cos (/.t + (3/2) tan-i(tr) - {^tT^R^ /2)/{l + fT^))]] . 

(C2) 

APPENDIX D: TIME EVOLUTION OF THE WAVE FUNCTION AND EXPONENTIAL 

DECAY OF SUPERPOSITIONS 

As mentioned in the introduction, an alternative form of the collapse equation has been given 
in (I), which differs from ()13p by a change of measure for the noise; see (35) and (37) of (I). The 
advantage of this alternative formulation is that it can be expressed in terms of a linear, but not 
norm-preserving, equation ( (34) of (I)), which is simpler to solve. Upon normalization and change 
of measure, one recovers the usual collapse dynamics. 

Let us specialize to the case of a single particle; let us moreover set // = 0, as we want to focus 
only on the collapse mechanics. Then, for the mass density coupling considered in this paper, the 
linear equation reads: 

d\x{t)) 



V7 / d^xM{x)(l){x,t) - 27 / d^x / d^yM{x)M{y)F{x - y,t) 



\x{t)) ■ (Dl) 



dt 

The random field (/)(x, t) is now supposed to be a Gaussian thermal field with respect to a new 
measure Q, having mean and correlator D{x — y,t — s). The relation between the statistical 
averages with respect to this measure and the averages with respect to the physical measure used 
throughout this paper (which we shall call P from now on) is: 

Ep[/(t)] = EQ[m{x{t)\xm , (D2) 

where f{t) is a generic random function of time. 
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Because of the special form ([3]) of the particle density operator M{x), (|Dip can be readily solved 
in the coordinate representation x(x, t) = {x\x{t)): 



x{^,t) = exp -^m<l>(x, t) — 2^171 1(0, t) x{^,0) 



(D3) 



with: 



$(x,i) 



ds(j)[x, s), 



1(0, t) 



dsF{6,s) 



(D4) 



'0 JO 

{I{x,t) has been first introduced in (|97p .) Let us fix an arbitrary time t. Then the random field 
^{x,t) is a Gaussian field in the variable x, with mean and correlator equal to: 



EQ[$(x,t)] =0, EqI^x, t)^y, t)] = 21 [x - y,t) 



(D5) 



The above statistical properties refer to the measure Q, while we need them to be expressed with 
respect to the physical measure P. ()D2p allows us to switch between the two measures, once the 
squared norm {x{t)\xi't)) has been computed. 

In analogy with the discussion of Sec. 3, let us consider an initial state of the form: 



Xix, 0) = ai6^{x - f ^)^/2 + a2d^{x - f^f^^ , 



(D6) 



corresponding to the superposition of two states well localized around r ^ and r ^ respectively. By 
substituting it into ()D3P and normalizing the wave function, one obtains for the collapse probabil- 
ities: 



pi{t) = |ai(t)|2 
P2{t) = |ai(t)|2 






Pi 



g2v/7m*{r i,i) _^ 2y7m'I>(r 2,t) 



(D7) 



with pi = I ail and p2 = \ct2\ ■ Using ()D2p . together with the equation 



{X{t)\xit)) = Pi exp 2^m$(f\t)-47m2/(0,t) 



+ p2exp 



2^m^{f^,t) - 47m^/(0, t) 



(D8) 



we can compute the average of the product pi{t)p2{t). 

Due to the statistical properties (IDSp . the joint probability density of the two random variables 
^{r^,t) and <l>(r^,t) reads: 



12 n / 2 12 



27rv^ 



exp 



Of (^(f ^ t))^ - 2atbt^{r',t)^{r^,t) + at(^(f ^ t)) 
2(a? - 6?) 



(D9) 
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with at = 21(0, t) and h = 21{f^ - f'^,t). Using now ([D2]), (JDT]), ^Bl and (JD9]) we get: 

x^ +y^ —2rtxy 



1 f+oo exp 

EF\piit)p2it)] = piP2e-^^'^ <"' == / dxdy 



8'ym'^at{l—r^) 



jz:;:T)+x + y 



STTjm'^aty^l -r1 J-oo pie^ + p2ey 

(DIO) 

with Tt = bt/at; we have also relabeled x = 2^m^{f^,t) and y = 2^m^{r^,t). ()D10|1 can be 

further simplified by making the change of variables t = {x + y)/2, s = {x — y)/2. In this case, the 

two integrals decouple and one gets: 

where T{t) = ^m'^at{l — Vf) corresponds to the definition (J19p . The final integral gives a finite 
contribution as T{t) -^ oo, which proves that the decay of the superposition is exponential in time, 



and proportional to e~^(*V\AX^- I^ particular, by using the inequality 

Pie''+p2e~'' > m(e'*+e~'') = 2mcoshs , (D12) 

with m = min{pi,p2} (here we assume that m 7^ 0; the trivial case m = can be treated 
separately), one has: 

+00 g-s2/4r(t) I r+oo ^ ^ 

ds < -r= / ds — -— = -r= ; (D13) 

-00 pie'^ + P2e ^ 2m J _^ coshs 2m 

Collecting all results, we can write: 

Ep[pi(i)P2(t)] < Ep[pi(0)p2(0)]-^^e-rW . (D14) 

APPENDIX E: UNPARTICLE THERMAL CORRELATION FUNCTIONS 

We take the unparticle thermal correlation function to be given by an average over thermal 
correlation functions for particles of mass /i > 0, using the same weighting function p(/x^) that is 

n 

used [25|] to generate the unparticle propagator from the propagator for a boson of mass fx, 

Writing the left hand side of (j94p as D{x,t,fj,) so as to explicitly show the mass dependence, the 
thermal unparticle correlation function Du is then given by 

/•oo 

Du{x,t)= d^?p{^?)D{x,t,^i) . (E2) 

JO 
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Substituting ([MD and ((95]), we thus get 

'■°° '■ (fk 1 



Du{x,t) =(d-l)A2(i-'^) / d/i2(^2)d-2 / ^LJ __ cos(fe-x)cos(a.fct) 

Jo J V^'^J '^fc g T — 1 

=(d - 1)A2(i-'^) I ^ cos (fc • x) j^ d^,\^?r'^^-^ cos {u,t) , 

(E3) 



where in the second hne we have isolated those factors of the integrand that exphcitly depend on 
^. Since oj'j, = k"^ + p?' ^ we can change integration variable in the inner integral from /i^ to ujf,, by 



using 

dl? = 2u;udi^k , (/x2)'^-2 = (.;2 _ ^2)d-2 ^ (E4) 

which gives 

/ d^2(^2)d-2 cos(u;fct)=2/ du; ^ ^_^ ^ cos (u;t) , (E5) 

Jo '^fe e~T - 1 JA: e^r- - 1 



where we have relabeled the dummy integration variable iOk as lo. Substituting this into (jESj) we 
get 

-^cos(A:-f) / duj ^\J ' cos(a;t) . (E6) 

v^^) Jk e T — \ 



Corresponding to this formula, the function Fij{x,t) introduced in (jll|) is given by 

'* ' d^^ _/r. ^\ f^ du {lo^ - k'^f-^ 

Ik ^ e^ — 1 



Fuix,t)= f dsDu{x,t-s) = 2{d-l)A'('~') fAh^cos{k-x) H ±^^^^-±1— sin (cot) , 
Jo J (2vr)^ Jk ^ e T - 1 

(E7) 

and the integral appearing in the rate function T{t) of (J17p is given by 

Iuix,t)^ ['dsFu{x,s)=2{d-l)A'('-''^ [-^cos{k-x) T ^ii^^!^Z^!t! [l -cos (^t)] . 
Jo J (2vr)^ Jk ^ e T - 1 

(E8) 

From (|E7p we can read off the Fourier transform defined in (j70p . from which the energy production 

is calculated through ([7T]) . 

roo T / 2 l.2\d-2 

4/(fe,t)=2(d-l)A2(i-^)/ ^ (^ ^:; ) sin (^t) . (E9) 

Jk ^ e^r- — 1 

Note that in all of these formulas, the scale parameter A appears as an overall factor, which then 
combines with the noise coupling 7 to give a new effective coupling 7A^''^~'^''. 

The correlation function Du{x^t) can be written in several alternative forms. Performing the 
angular average over fc, we get 

Duix, t) = [d- 1)A2(i-^) / ^ sin (A;|x|) / dcu ^^_^ ^ cos (cut) , (ElO) 

JO ^ pl Jk e'^r- — 1 
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which on interchange of orders of the k and uj integrations becomes 

Jo e~^r- — 1 JO ■^ fI 

Making the change of integration variable k = cov, this can be further rewritten as 



u2\d-2 



a^d) r . ^2(^^-1) COS (tjt) [^ dvv . , _,.,^ 2^d-2 

^ ' \ du) j3^ !^ — ^ / -^—r sm {vuj\x\) {1 - v'^y ^ 



(Ell) 



Z5^(x,t) = (d-l)A'^i-^M do; —^ ^^ / ^:^sin(«a;|f|)(l-i>^)'*-' . (E12) 

-'0 e~^r- — 1 Jo ^ fI 

The integral over i; in ()E12p converges only for Reti > 1. However, by an integration by parts this 

integral is transformed as follows, 

^' ^ sm (.elf I) (1 - v-r- = £ ^J^^ cos (..|x I) (1 - vr-' , (E13) 

which gives an analytic continuation around the simple pole at d = 1, expressed in terms of a f 
integral that now converges for Red > 0. Substituting (1E13|1 into (|E12p gives a formula for the 
correlation function which is now manifestly finite for Re d > 0, 

Du{x,t) = 1a^(^-^) rd ^Jr^"'^ [''-^,cos{vu\x\)il-vr-' . (E14) 

^ JO e^r- — 1 JO "^ 

The corresponding formulas for Fu{x, t) and Iu{x, t) are now obtained by the replacement of cos(a;i) 

by sm(Lot)/Lj and [1 — cos{uit)]/uj'^, respectively, 

Fu{x,t) = iA^(i-^) rd ^Jr^""'^ f%.os{vu:\x\){l-vr-' , (E15) 

^ JO e~^r~ — 1 JO ■^ 

and 

Iu{x,t) = lA^a-'^) r, y-'[^--o<-t)] ["dVs{vu\mi-vy-' . (E16) 

2 Jo e T - 1 Jo -^ 

From ()E16P we find for the subtracted integral that enters into T{t), 

^ JO e^^ — 1 JO ^ 

(E17) 

giving an expression that is manifestly positive. 

Let us now return to (|E9p . and use it to calculate the energy production. Substituting (|E9p 

into (j7ip we get, for a single particle with mass-coupled unparticle noise, 

d^,, ,, 7mA2(i-'^)(d-l) f^ „ , f^ dio iuo'^ - k'^f-^ , . ,^ , 

-TvHp{t) = l ^ '- dkk^ — ^ ._c ' sin(a;t) , (E18) 

at TT Jo Jk ^ e~^r~ — 1 



which on interchanging the orders of the k and a; integrations becomes 

,, 7mA2(i-'^)(d-l) /-"^ dw sin (cjt) /"'^ , ,, ., ..^^ ., 



Axri/oft) = 7^A^(^-'^n^ - 1) r d^^in^ , ^,^ , ,^ ^ .2^d-2 
dt ^^' 7r2 J^ uj ^^ _i Jo 



42 

Making the change of variable k = cou^'"^, dk = {\/2)Lou~^''^du in the inner integral, it can be 
evaluated in terms of the Euler B function; then using [d— l)T{d— 1) = T{d) one gets the compact 
expression 

dt '^^ (2i)2 r(3/2 + (!) y„ eT* - 1 

with the integral convergent for Red > 0. Integrating over t, one gets the corresponding formula 
for the total energy production, 

^ .. . N ^ rr , . 37mA2(i-««) r(3/2)r(d) f^ , a;^^"^ [l - cos (wt)] ,^^^, 
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